Astronomy & Astrophysics manuscript no. 15938 


©ESO 2010 


December 6, 2010 





The Westerbork Hydrogen Accretion in Local Galaxies 

(HALOGAS) Survey 

I. Survey Description and Pilot Observations 

George Heald\ Gyula J6zsa\ Paolo Serra^, Laura Zschaechner^, Richard Rand^, Filippo Fraternali^, Tom Oosterloo^'"^, 

Rene Walterbos^, Eva Jiitte^, and Gianfranco Gentile' 



o 

<N 
O 

Q 



O 
U 

O 

03 



Netherlands Institute for Radio Astronomy (ASTRON), Postbus 2, 7990 AA Dwingeloo, the Netherlands 

University of New Mexico, 800 Yale Blvd, Albuquerque, NM, USA 

Astronomy Department, Bologna University, Via Ranzani 1, 40127 Bologna, Italy 

Kapteyn Astronomical Institute, Postbus 800, 9700 AV Groningen, the Netherlands 

Department of Astronomy, New Mexico State University, RO. Box 30001, MSG 4500, Las Gruces, NM 88003, USA 
Astronomisches Institut der Ruhr-Universitat Bochum, Universitatsstr 150, 44780 Bochum, Germany 
Sterrenkundig Observatorium, Ghent University, Krijgslaan 281, S9, B-9000 Ghent, Belgium 



Received 15 October 2010; accepted 27 November 2010 



ABSTRACT 



We introduce a new, very deep neutral hydrogen (Hi) survey being performed with the Westerbork Synthesis Radio Telescope 
(WSRT). The Westerbork Hydrogen Accretion in LOcal GAlaxieS (HALOGAS) Survey is producing an archive of some of the 
most sensitive H i observations available, on the angular scales which are most useful for studying faint, diffuse gas in and around 
nearby galaxies. The survey data are being used to perform careful modeling of the galaxies, characterizing their gas content, mor- 
phology, and kinematics, with the primary goal of revealing the global characteristics of cold gas accretion onto spiral galaxies in the 
local Universe. In this paper, we describe the survey sample selection, the data acquisition, reduction, and analysis, and present the 
data products obtained during our pilot program, which consists of UGC 2082, NGC 672, NGC 925, and NGC 4565. The observations 
reveal a first glimpse of the picture that the full HALOGAS project aims to illuminate: the properties of accreting H i in diff'erent types 
of spirals, and across a range of galactic environments. None of the pilot survey galaxies hosts an Hi halo of the scale of NGC 891, 
but all show varying indications of halo gas features. We compare the properties of detected features in the pilot survey galaxies with 
their global characteristics, and discuss similarities and diff'erences with NGC 891 and NGC 2403. 
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It is becoming increasingly clear that gas accretion is a crucial 
piece of the puzzle of galaxy evolution. Accretion has been in- 
voked as a possible solution to several problems in galactic as- 
trophysics. A typical accretion rate of 1 M© yr"^ seems to be 
necessary in order for spiral galaxies to replenish gas consumed 
in the star formation process - it has been noted that otherwise, 
galaxies w ould use up the ir gas on timescales shorter than their 
ages (e.g., Larson et al.l 1980J. The constancy of the SFR in 
the solar neighborhood ( |Twarog||1980[ Binney et al. |200Q| ) ar- 
gues that this inflow of fresh gas should locally be fairly con- 
stant. Similar accretion rates of metal-free gas m ay be needed 
to match the chemical evolution of the disk (e.g., (Lacey & Fall] 



1985^ 



Schonrich & Binney 2009). Accretion of fresh gas has 
also been invoked to explain the "G-dwarf prob lem," i.e. th e 
small number of local stars with low metallicity ('Pagel""1997 ), 
for instance by Chi appini et al.| ( pW7] ). Large scale disturbances 
in galaxies, such as lopsidedness (both morphological and kine- 
matical) and warps, may be the result of accretion ( |Boumaud| 
let al.||2005| [Sancisi et aT]|2008| ). In their recent review of the 
observational evidence for cold accretion, Sanci si et al.| ( |2008| ) 
infer a global average observed neutral gas accretion rate onto 



galaxies of > 0.2 M© yr"^ Such accretion may take place in two 
ways: through the merging of gas-rich satellites, and "pure" gas 
accretion from the intergalactic medium (IGM). 

On theoretical grounds, a mode of accretion referred to as 
"cold" should continue in galaxies in low density environments 
through z = (e.g., Keres et al. 2005 ). Note that this usage of 
the word "cold" diff'ers from the meaning typically assumed by 
observers: in the theoretical picture of accretion, cold gas is that 
which is never shock heated to the virial temperature of the halo 
as it falls into the galactic potential, reaching maximum tempera- 
tures below 10^ K during the accretion process; observationally, 
cold gas means gas at temperatures up to nearly 10"^ K. Despite 
this diff'erence in nomenclature, it seems clear that a significant 
portion of the accreting gas in local galaxies is expected to be at 
relatively low temperature, and that this gas phase will be traced, 
at least in part, by H i. 

Indeed, deep observations of a few nearby spirals - most 
notably the edge-on NGC 891 ( Qoste rloo^et al.||2Q07| ), 



the 



moderately-inclined NGC 2403 (,Fratemali et al.|2QQ2| ),"and the 
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face-on NGC 6946 ( |Boomsma et al.||2Q08| ) - have shown that 
large quantities of cold gas exist in their halos, at large vertical 
distances from the thin disk. In the case of NGC 891, observa- 
tions with the Westerbork Synthesis Radio Telescope (WS RT) 
revealed a huge neutral gas halo (e.g. Qosterloo et al.|2Q07 see 
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their Figure 1) which extends up to a vertical distance of 22 kpc, 
with massive filamentary structures far from the disk. In total, 
the Hi halo of NGC 891 contains almost 30% of the galaxy's 
total Hi mass. 

How can so much gas be located at such large vertical dis- 
tances from the disk? Star formation-driven disk-halo flows, as 
described by galactic fountain (e.g., Bregman 1980) or chimney 
( Norman & Ikeuchi,1989j models, are likely part of the explana- 
tion. The morphology of the extraplanar gas seems to be related 
to the properties of the underlying star forming disk (at least in 
its ionized phase; e.g.. Rand 1996 ), and whether ionized gas is 
found in the halo or not seems to depend on the energy density 
provided by star formation ( Rossa & Dettmar|2003 ). 

On the other hand, it has been shown that the rotational ve- 
locity of the halo gas (both in its neutral and ionized phases) 
decreases linearly with increasing distance from the midplane 
( [Fratemali et"arf 2005 : Hea ld et al. 12006 : Kamphu Tset al.|2007| ). 
This kinematic "lag" - which allows us to distinguish halo gas 
from disk gas in moderately inclined galaxies such as NGC 2403 
- is much steeper than can be accounted for with simple galac- 
tic fountain mod els (the d iscrepancy is observed to be a factor 
of two or more; Heald et al. 2007 ). A good deal of theoretical 
efl'ort has been put into this problem, and a range of more so- 
phisticated models to explain the kinematics of gaseous halos 
have been put forward (using A/^-body-F smoothed particle hydro- 
dynami c (SPH) s imulations , 'K aufmann etal." f20 06j); hydrostatic 
models, 'Benjamin ( 2000 ); Ba mabe et al. (2006 ) ; Marinacci et al. 

tlOb); SPH including spiral arms and bars, [Struck & Smith 
)9 ); and cloud-corona interaction, [ Marinacci et al.| ( |201 0a')). 
Such efl'orts have shown how the steep lag in NGC 891 might 
occur. Additionally, a constant accretion of low-angular momen- 
tum gas has been invoked, in combination with the simple galac- 
tic fountain picture, to slow the halo gas (the required accretion 
rate would be 3 M© yr"^ in NGC 891 based on modeling work 
by |Fraternali & Binney ( 2008| ) - in this picture, 90% of the halo 
gas originates from galactic fountain, and 10% from accretion). 
Distinguishing between these possible models requires a larger 
sample of gaseous halos and measurements of their kinematics. 

It has been pointed out that if the largest H i filament in the 
halo of NGC 891 was created by disk star formation, an energy 
equivalent to ~ 10^ supernova e would be required. T his sug- 
gests an external origin instead ( Qosterloo et al.||2007| ). A final 
piece of evidence suggestive of accretion in NGC 891 comes in 
the form of small, 10^ M© clouds which are counter-rotating in 
the halo. A filament wit h a mass of 10^ M© is found in NGC 2403 
(Fratemali et al. ]|2002| , and a counter-rotating -10^ M© cloud 



is a lso detected in the extraplanar region of NGC 5746 ( |Rand & 
[Benjamin 2008) . 

Tracers of ongoing accretion can also be found closer to 
home: many of the MW's high velocity clouds (HVCs) have 
low metallicity ( [Tripp etaL 20031; together wi th the increas 



ingly tight constraints on their distances ( Wakker et al. 2007 



|2008) , the evidence is consistent with the idea that at least part 
of the HVC population is inf ailing primordial gas. HVC analogs 
have also been detected in the outskirts of M3 1 ( [Thilker et al.| 
[2004| ). In both the MW and M31, stella r streams are observed 
via de ep optical observations (Belokur ov et al.[[2007[ Ferguson 
2007) ). The H i filament in the halo of NGC 891 may be the cold 



gas counterpart of this type of stream. It has been recently rec- 
ognized that NGC 891 also has a giant stellar stream looping 
around the disk, pointing to recent interaction (Mouhcine et al. 
2010| ). Observing such features in H i emission has the advantage 



Halo gas is also present in at least one low surface bright- 
ness (LSB) galaxy (NGC 4395) with very little ongoing star for- 
mation, and a 10^ M© gas cloud in that system is l ikely in the 
process of being accreted ( Heald & Oosterloo|2008 ). Moreover, 
evidence suggests the presence of a neutral gas halo in the su- 
perthin LSB UGC 7321 ( [Matthews & Wood|2003| ). On the other 
hand, there is no evidence for large amounts of H i in the halo of 
the very massive NGC 5746 ( Ra nd & Benjamin|2008| ). 

Why do some galaxies have massive gaseous halos, but not 
others? Is the Hi halo of NGC 891 typical or exceptional? At 
present, the extremely limited number of sufficiently deep ob- 
servations restricts our ability to draw general conclusions about 
the prevalence of accretion in the local Universe. A compre- 
hensive investigation is needed, using a representative sample 
of spiral galaxies. Most of the existing deep observations were 
prompted by earlier shallow observations which already showed 
evidence for an exceptional halo population or other peculiari- 
ties. Moreover, the targeted galaxies were chosen individually, 
rather than as a part of a larger program to sample a wide range 
of galactic disk properties such as star formation rates, warps, 
and lopsidedness. 

As a first attempt to provide such information, we have 
started the Westerbork Hydrogen Accretion in LOcal GAlaxieS 
(HALOGAS) Survey. With this project we have the goal to make 
a significant step forward in understanding the process of cold 
gas accretion in spiral galaxies, by investigating the general phe- 
nomenon for the first time. 

The HALOGAS Survey observations are designed to detect 
and characterize neutral extraplanar gas by measuring its total 
mass, rotation speed (in edge-ons), radial inflow or outflo w (in 
inclined ga laxies such as NGC 2403; see the work by Frate mali] 
et al.|2002 ), and distribution (smooth, filamentary, or clumpy?). 
Deriving the extraplanar gas kinematics is of crucial importance 
- these parameters are key predictions of halo gas and accretion 
models. Comparisons between the extraplanar gas characteris- 
tics and the properties of the spiral disk should clarify the roles 
of star formation and accretion. 

The HALOGAS data will form a publicly available repos- 
itory of the deepest Hi observations of nearby spiral galaxies. 
Additional studies which will be enabled by this dataset are nu- 
merous, and include: 

1 . Search for background continuum source variability at low 
(~ 1 mJy) flux levels, over the full HALOGAS survey area 
(totalling ~ 6.2 square degrees within the WSRT primary 
beam FWHM); 

2. Detailed analyses of the gravitational potential by measuring 
rotation curves extending far into the faint outer disks, and 
additionally by using the combination of the gas dispersion 
and vertical scaleheight; 

3. Studies of warps and lopsidedness in spiral galaxies; and 

4. The details of the Hi cutoff at the edge of the gaseous disk, 
at sens itivity levels comparable to the requirement specified 
by e.g. [Corbelli & Bandiera[ ( [20Q2| ) to examine the H i - H ii 
transition. 

The survey observations are being completed at the WSRT. 
Our first semester of observations was designed as a small pilot 
program, and that pilot survey is the focus of this paper. 

The paper is organized as follows. We describe the selection 
of the survey sample in §|2] The observational setup and data 



reduction path are outlined in 
Pilot Survey are presented in §[4 



that we directly obtain kinematical information. 



3] Results from the HALOGAS 
We discuss the first results from 



the HALOGAS Survey in §|5] Finally, multi- wavelength supple 
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ments to the HALOGAS Survey, and our data release plans, are 
discussed in §[6]and[7l 



2. The HALOGAS Sample 

To reach the sensitivity required to detect and characterize faint, 
diffuse H i, with characteristic column densities of a few times 
10^^ cm"^, each of the HALOGAS targets is being observed at 
the WSRT for a total of 10 x 12hr. We would like to have an un- 
biased look at the accretion process in the local Universe. From 
this perspective, non-detections of extraplanar gas are as inter- 
esting as the detections - by restricting attention to galaxies al- 
ready suspected to show clear signs of accretion, we would learn 
less about the ubiquity of the accretion process in the spiral pop- 
ulation, or the relationship between galaxy parameters and the 
accretion process. We therefore deliberately avoided selecting or 
rejecting sample galaxies based on existing Hi observations, or 
other subjective judgements about how interesting a target might 
appear when considered as an individual target. The final sam- 
ple of galaxies was, ultimately, chosen on the basis of a small 
number of basic and neutral selection criteria, each of which are 
designed to maximize our ability to detect and study the prop- 
erties of any diff'use gas which might be present in the target 
systems. 

To select the galaxies for our sample, we searched in the 
Tully|(p^88]) catalog, requiring the following characteristics: 



(i) The sample is made up of barred and unbarred spirals with 
Hubble types between Sa and Sd. 

(ii) Because the WSRT is an east- west interferometer, we require 
that the source declination S > -\-25° to ensure that the syn- 
thesized beamsize is < 30''. 

(iii) We only select nearby galaxies with distance D < II Mpc, 
in order to observe at adequate physical resolution to resolve 
the morphology of the galaxies (for instance, the character- 
istic WSRT beamsize of 15'' = 0.7 kpc at a target distance of 
10 Mpc). 

(iv) In order to ensure that each galaxy is sampled by a reason- 
ably large number of angular resolution elements, we re- 
quire a large optical isophotal diameter at 25th magnitude 
D25 > 3'. Typically, Hi disks are larger than D25. 

(v) Finally, we require that the systemic velocity is > 100km s"\ 
in order to avoid confusion with MW H i. 



Within these criteria, we selected all intermediately inclined 
(50° < / < 75°) and edge-on (/ > 85°) galaxies. Galaxies in the 
former category reveal connections between halo gas and star 
formation, and allow investigation of rotation and radial motions 
in the halo gas (e.g. Fr aternali et al. 2002) ; edge-on galaxies are 
useful for learning about the extent of halos, and the vertical 
change in rotation speed of the halo gas (e.g. |Oosterloo et al. 
2007| ). Visual inspection was utilized only to discard galaxies 



with incorrectly catalogued inclination angles. We also excluded 
from our new observations two galaxies which meet our criteria, 
but which are already deeply observed in Hi (namely, NGC 891 
and NGC 2403). We will incorporate the properties of these 
galaxies into our final statistical sample since they match the 
sample criteria described above. The final sample of 24 galaxies 
is Hsted in Table[T] The columns are (1) the UGC ID; (2) NGC 
or Messier IDs if available; (3) the Hubble type; (4) the distance 
Jtss in Mpc from T ully| ( p^88 ); (5) the adopted distance Jbest 
in Mpc; (6) the method for determining Jbest; and our best esti- 
mate using the tabulated method; (7) the systemic velocity v. 



the blue magnitude Mg; (11) the rotational velocity Vrot; (12) the 
star formation rate (SFR); and (13) references used in deriving 
the SFR values in column (12), see below for details. The UGC 
IDs of the four targets which have been observed in the pilot 
survey described in this paper are indicated with a bold font. 
The galaxy parameters listed in Table [T] are tak en from|Tully 



(1988), except for Vrot, which is from HyperLEDA ( [Paturel etaL 



2003 ), the adopted distance Jbest. and the SFR, both of which are 
discussed below. Given the broad range of Hubble type, SFR, 
mass, and environmental parameters (see Table [2] and the dis- 
cussion below) exhibited by nearby spiral galaxies, this sample 
size stikes an acceptable balance between the number of targets 
required to reasonably sample the parameter space, while still 
resulting in a total amount of observing time which does not be- 
come prohibitive. 

One of the sample selection criteria is to select only galax- 
ies closer than 11 Mpc. This selection is based on galaxy dis- 
tances taken from |Tully| ( |1988 ), which we indicate as Jtss- More 
recent distance estimates are available for all selected galax- 
ies. Adopting the best one is important when measuring, e.g., 
mass and linear size of the detected Hi structures. Therefore, 
we searched the NASA/IPAC Extragalactic Database (NEdQ 
distance database to obtain the most accurate distance of each 
galaxy, which we indicate as Jbest- Table [T] lists both Jt88 and 
Jbest for all galaxies in the sample. We define Jbest using one of 
the following methods, sorted according to decreasing priority, 
depending on which measurements are available: 

1 . Median of all Cepheid and/or tip of the red-giant branch dis- 
tances: used for 10 galaxies. 

2. Median of all distance estimates obtained with any of the 
following methods: planetary-nebulae luminosity function, 
brightest stars, radio and optical SNII, masers, globular- 
cluster luminosity function: used for 5 galaxies. 

3. Median of all Tully-Fisher distances (including Jtss)- used 
for 8 galaxies. 

4. ' ~ 



TuUy ( 1988) estimate: used for 1 galaxy. 



It is worth noting that Jbest is significantly diff'erent from Jt88 
(> 2 Mpc) in about a third of all galaxies. For NGC 3198 and 
NGC 4414 Jbest is obtained with method 1, and is respectively 
3.7 and 8.1 Mpc larger than Jtss- For NGC 4062, NGC 4274, 
UGC 2082, UGC 4278 and UGC 7774, 4est is obtained with 
method 3 and is always larger than Jt88 (the last galaxy is the 
most extreme case, with a diff'erence of 17.6 Mpc between the 
two estimates). In general, Jtss is always one of the lowest 
(when not the lowest) distance estimates available. 

The method that we use to calculate star formation rates 
also bears some discussion. The galaxies that make up the 
HALOGAS sample are not otherwise included in a uniform ob- 
servational sample at other wavelengths, in particular the in- 
frared. We therefore use diff'erent methods for individual tar- 
gets to calculate SFR, depending on what data are available. We 
pref er to use the co mbination of Ha and TIR luminosity given 
by Kennicutt et al.| ( [200 9), where Ltir is either calcula ted from 
IRAS or Spitzer fluxes via the prescriptions given by |Dale & 
Helou (2002'). The IRAS form is used where possible, and the 
Spitzer form where IRAS fluxes are not available or provide only 
upper limits for the 25/dm flux. IRAS fluxes are obtained from 
the references listed in Table [T] and the Spitzer fluxes are from 
the Local Volume Legacy (LVL;||Dale et al. 2009). Ha fluxes are 



sys 



^ The NASA/IPAC Extragalactic Database (NED) is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under 



in km s ; (8) the inclination /; (9) the optical diameter D25 ; (10) contract with the National Aeronautics and Space Administration. 
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Table 1. The Full HALOGAS Sample. 











^best 


^best 


^sys 


i 


D25 


Mb 


Vrot 


SFR 


SFR 


UGC 


Other IDs 


Type 


(Mpc) 


(Mpc) 


method 


(kms-1) 


n 


(arcmin) 


(mag) 


(kms-^ 


(Moyr-0 


refs 


1256 


NGC 0672 


SBcd 


7.5 


7.6 


2 


425 


70 


6.4 


-18.65 


130.7 


0.22 


HI, 12 


1831 


NGC 0891 


SAb 


9.6 


9.2 


1 


529 


84 


12.2 


-19.96 


212.2 


2.4 


HI, 14 


1913 


NGC 0925 


SABd 


9.4 


9.1 


1 


554 


54 


11.3 


-19.66 


102.4 


0.83 


HI, 14 


1983 


NGC 0949 


SAd 


10.3 


11.3 


2 


610 


52 


3.5 


-17.85 


90.9 


0.26 


H1,I1 


2082 


- 


SAc 


10.7 


14.4 


3 


710 


89 


5.8 


-18.55 


86.6 


0.023 


-,I1 


2137 


NGC 1003 


SAcd 


10.7 


11.6 


3 


626 


67 


6.3 


-18.61 


95.5 


0.34 


HI, 12 


3918 


NGC 2403 


SAcd 


4.2 


3.2 


1 


132 


62 


23.8 


-19.68 


121.9 


1.0 


HI, 14 


4278 


IC 2233 


SAd 


10.6 


13.6 


3 


565 


90 


4.3 


-17.45 


79.2 


0.11 


HI, 15 


4284 


NGC 2541 


SAcd 


10.6 


12.0 


1 


553 


67 


7.2 


-18.37 


92.1 


< 0.27^ 


H2,I2 


5572 


NGC 3198 


SBc 


10.8 


14.5 


1 


660 


71 


8.8 


-19.62 


148.2 


0.61 


H2,I4 


7045 


NGC 4062 


SAc 


9.7 


16.9 


3 


769 


68 


4.5 


-18.27 


140.5 


0.22 


H2,I2 


7322 


NGC 4244 


SAcd 


3.1 


4.4 


1 


247 


90 


15.8 


-17.60 


89.0 


0.058 


HI, 15 


7353 


NGC 4258 (M106) 


SABbc 


6.8 


7.6 


1 


449 


71 


17.1 


-20.59 


208.0 


1.4 


HI, 13 


7377 


NGC 4274 


SBab 


9.7 


19.4 


3 


922 


72 


6.5 


-19.22 


239.9 


0.31 


H3,I2 


7539 


NGC 4414 


SAc 


9.7 


17.8 


1 


720 


50 


4.5 


-19.12 


224.7 


1.3 


H2,I4 


7591^ 


NGC 4448 


SBab 


9.7 


9.7 


4 


693 


71 


3.8 


-18.43 


221.6 


0.056 


-,I2 


7766 


NGC 4559 


SABcd 


9.7 


7.9 


3 


816 


69 


11.3 


-20.07 


113.4 


1.1 


H4,I4 


7772 


NGC 4565 


SAb 


9.7 


10.8 


2 


1228 


90 


16.2 


-20.34 


244.9 


0.54 


H5,I4 






SAd 


6.8 


24.4 


3 


526 


90 


3.5 


-15.57 


79.4 


0.0074 


HI, 15 


7865 


NGC 4631 


SBd 


6.9 


7.6 


1 


613 


85 


14.7 


-20.12 


138.9 


2.1 


HI, 14 


8286 


NGC 5023 


SAc 


6.0 


6.6 


1 


400 


90 


6.8 


-17.29 


80.3 


0.032 


HI, 15 


8334 


NGC 5055 (M63) 


SAbc 


7.2 


8.5 


3 


497 


55 


13.0 


-20.14 


215.5 


1.5 


HI, 14 


8550 


NGC 5229 


SBc 


6.4 


5.1 


2 


365 


90 


3.5 


-15.82 


57.3 


0.014 


HI, 15 


9179 


NGC 5585 


SABd 


7.0 


8.7 


2 


303 


51 


5.5 


-17.96 


79.1 


0.26 


H1,I1 



^ The S FR value for NGC 2541 is strictly speaking an upper limit, because the IRAS ISfim flux is catalogued as a non-detection by lMoshir] 

[etar(1990). 
^ No H I observing time granted yet. 



References. (HI) Kenn icutt et al. 'f2008^; ( H2) [Moust akas & Kennicutt'f2006^; (H3)"Hameed & Devereux"f2005^; (H4) 'Ken nicutt et al.|([2009l); 
(H5) Robitaille et al.,(2007| ); (Il) |Lisenfeld et al.| ( |2007| ); (I2) .Moshir et al.. ( .1990j ; (I3) ,Rice et al.. ( .1988j ; (I4) .Sanders et al.. ( .2003| ); (I5) |Dale et al.] 
( [2009F : 



obtained first from 'Ken nicutt et aL ( 2QQ8| ), and for targets not 



included in that work, the literature was searched for fluxes. In 
all cases luminosities are calculated using the adopted distances 
listed in Table [T] In two cases, UGC 2082 and NGC 4448, no 
Ha fluxes are available; but since IRAS 60yum and lOOyum fluxes 
are available, we calculate Lfir using the equation given by 
Rice et al. ( 1988 1, and use SFR(Lfir) as described byl Kenn icuttj 
( |1998| ). The references given in Table [T] were used to obtain the 
Hof and IR fluxes. The derivations of the SFR for the HALOGAS 
sample will become more uniform with the release of the Wide- 
field Infrared Survey Explorer (WISE; Wrightj2008j) data, which 
will also be well matched in angular resolution to our WSRT ob- 
servations. 

In Table [2j we show environmental information about the 
target galaxies. The data are from the group cata log developed 
from the Nearby Optical Galaxies (NOG) sample ( [Giuricin et al. 
|200 0). The columns in the table are: (1) galaxy ID; (2) number 
of group members (using the percolation- 1 grouping algorithm, 
see below); (3) NOG group (NOGG) number; (4) total group B 
luminosity; (5) distance from group center of mass to outermost 
member; (6,7) minimum and maximum group member velocity 
relative to group center of mass; (8) contribution of the target to 
the group B luminosity; (9) fraction of maximum group member 
B luminosity; (10) normalized distance to group center of mass; 
(11) target velocity relative to group center of mass. 

We have chosen to use group definitions resulting from the 
Giuricin et al.| ( |2000| ) "PI" algorithm (which is a percolation 
friends-of-friends algorithm that keeps the velocity grouping 



condition fixed irrespective of the galaxies' redshift). Because 
this definition of groups is based on angular and velocity coinci- 
dences, there will be spurious groupings or separations. We note 
corrections to the group definitions where appropriate in the text, 
and otherwise use the information as a means to illustrate the 
overall environmental characteristics of the HALOGAS sample. 
The survey galaxies sample a fairly broad range in environment. 
One quarter (6/24) of the HALOGAS sample reside in the same 
large group (Coma I). Three targets are members of pairs. In 
only one case is there a galaxy with a strongly interacting neigh- 
bor: NGC 672, which was included in the pilot survey described 
in this paper. The range of galaxy properties with respect to the 
group to which they belong also exhibits a wide range of val- 
ues. The HALOGAS sample contains galaxies which dominate 
their group, as well as those which are minor group members 
and companions. If these smaller targets are gas rich, we may be 
probing gas accretion donors as well as recipients. Considering 
all of these properties, after observations of the full HALOGAS 
sample, we expect to be able to assess whether there is an envi- 
ronmental dependence on target galaxy halo properties. 

Because the HALOGAS galaxies are all nearby, there is the 
added advantage that most have been extensively observed at 
other wavelengths, for example as part of recent surveys such as 
the Spitzer Infrared Nearby Galaxies Survey (SINGS; Kennicutt] 
et al. 200 31) and the 1 1 Mp c Ha and Ultraviolet Galaxy Survey 
(1 lHUGS; ]Lee et alT 2007 ). The star formation process is clearly 
of interest with relation to the extraplanar, halo, and accreting 
gas that we are tracing, and we will therefore benefit by having 
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Table 2. Environmental properties of the HALOGAS sample galaxies. 



Galaxy 


# group 


NOGG 


-^5,group 


-^group 


^min,group 


^max,group 


L'B 1 ^5,group 


L'bI Lb, max 


R/ ^group 


V 


ID 


members 


number 


[10^0 Mo] 


[Mpc] 


[km/s] 


[km/s] 








[km/s] 


NGC 0672 


2 


88 


1.3 


0.01 


-52 


29 


0.64 


1.00 


0.55 


29 


NGC 089 P 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


NGC 0925 


6 


149 


3.1 


0.56 


-16 


47 


0.69 


1.00 


0.25 


-16 


NGC 0949 


6 


149 


3.1 


0.56 


-16 


47 


0.15 


0.21 


1.00 


47 


UGC 2082^ 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


NGC 1003^ 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


NGC 2403 


2 


319 


1.8 


0.24 


-14 


1 


0.88 


1.00 


0.13 


1 


UGC 4278 


6 


345 


1.4 


0.77 


-90 


35 


0.30 


1.00 


0.44 


21 


NGC 2541 


6 


345 


1.4 


0.77 


-90 


35 


0.24 


0.79 


0.21 


35 


NGC 3198 


2 


462 


1.7 


0.43 


-25 


1 


0.96 


1.00 


0.05 


1 


NGC 4062 


27 


63 r 


21.0 


1.27 


-163 


525 


0.03 


0.07 


0.86 


45 


NGC 4244 


16 


636 


3.1 


0.80 


-93 


83 


0.22 


0.66 


0.31 


-42 


NGC 4258 


21 


644 


9.0 


1.07 


-103 


160 


0.51 


1.00 


0.24 


-40 


NGC 4274 


27 


63 r 


21.0 


1.27 


-163 


525 


0.04 


0.10 


0.45 


194 


NGC 4414 


27 


63 r 


21.0 


1.27 


-163 


525 


0.05 


0.13 


0.21 


-3 


NGC 4448 


27 


63 r 


21.0 


1.27 


-163 


525 


0.02 


0.06 


0.37 


-79 


NGC 4559 


27 


63 r 


21.0 


1.27 


-163 


525 


0.12 


0.30 


0.42 


79 


NGC 4565 


3 


648 


14.9 


0.26 


-18 


107 


0.78 


1.00 


0.25 


-18 


UGC 7774 


21 


644 


9.0 


1.07 


-103 


160 


0.01 


0.01 


0.74 


3 


NGC 4631 


27 


63 r 


21.0 


1.27 


-163 


525 


0.40 


1.00 


0.30 


-91 


NGC 5023 


7 


723 


8.0 


0.45 


-94 


131 


0.04 


0.09 


0.68 


-77 


NGC 5055 


7 


723 


8.0 


0.45 


-94 


131 


0.37 


0.79 


1.00 


8 


NGC 5229 


7 


723 


8.0 


0.45 


-94 


131 


0.01 


0.02 


0.97 


-94 


NGC 5585 


8 


766 


3.1 


0.54 


107 


75 


0.07 


0.09 


0.52 


75 



^ Not included in the|Gi uricin et al.||2000| ) catalog due to low galactic latitude. 
^ No environmental information available; UGC 2082 is isolated. 
^ NOGG number 63 1 corresponds to Coma I. 



these data available for our investigation. Although the available 
data will help to characterize disk star formation rates, we note 
that the coverage is not enough to cover the full WSRT primary 
beam. Furthermore, we are pursuing deep observations at other 
observational facilities to fill gaps where necessary (see §[6]). 

3. Observations and data reduction 

The data are being obtained using the WSRT, which is an east- 
west interferometer consisting of 10 fixed (of which one is not 
used) and 4 movable 25 -m antennas. Baselines range from 36 m 
to 2.7 km, and the fixed antennas are on a regular grid with a 
spacing of 144 m. We used the array in its IMaxishort configura- 
tion, which is designed to optimize the imaging performance for 
extended sources in individual tracks. The correlator was set up 
to provide two linear polarizations in 1024 channels over a 10 
MHz bandwidth centered at the systemic velocity of each target. 
Each galaxy was observed for 10 full 12-hour tracks, and each 
track was bracketed by standard flux calibrator sources. 

Data reduction was performed using usual techniques in the 
Miriad software package ^Sault et al.|1995| ). A standard data re- 
duction path was used for most targets; exceptions were nec- 
essary due to individual circumstances and will be described 
where appropriate. First, the raw data were visually inspected, 
and samples corrupted by radio frequency interference (RFI) 
were flagged where necessary. System temperature variations, 
which are measured and recorded at each antenna during the 
observations, were corrected for in both the calibrator and tar- 
get data sets. Next, the standard calibrator sources were used to 
set the flux scale, as well as define and correct for the band- 
pass response. The radio continuum emission in the field was 
subtracted in the visibility domain and saved as "channel-0" 



data. An imaging and self-calibration loop was used on the 
continuum data to solve for time- variable antenna gain phases. 
The resulting antenna gains were subsequently applied to the 
continuum-subtracted spectral line data. The H i data were im- 
aged with a variety of weighting schemes (described below), 
creating several diff'erent data cubes for each target (see Table[3]). 
Offline Hanning smoothing led to a final velocity resolution of 
4.12kms"\ although we may choose to degrade this resolution 
in post-processing in order to maximize sensitivity to features 
with broad velocity width. Finally, Clark CLEAN deconvolution 
was performed in two stages, the first within masked regions 
containing the brightest line emission, and the second within the 
entire imaged region. 

Diff'erent weighting schemes are used in inverting the data 
from the u, v plane into the image domain. In practice, within 
miriad's invert task, we use a robust parameter of for in- 
termediate resolution and sensitivity. To maximize sensitivity to 
faint extended emission we additionally use a Gaussian u, v ta- 
per corresponding to 3C in the image plane. See Figures [T] [3] [5] 
andd 

For the pilot survey galaxies presented in this paper, we list 
in Table [3] a summary of the data cubes which have been pro- 
duced and analyzed for the purposes of the presentation in §[4] 
The columns are (1) Galaxy ID; (2) cube description; (3,4) beam 
size in arcsec and in kpc (note, the beam position angle is in all 
cases very close to 0° because WSRT is an east- west array); (5) 
rms noise level; (6) the 3cr column density sensitivity to emis- 
sion with a velocity width of 16 km s"^ ; and (7) the H i mass for 
a 3cr detection of a spatially unresolved source with the same 
16kms"^ velocity width. The Hi masses are calculated using 
Jbest from Table [T] 
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Table 3. Properties of Pilot Survey Datacubes 



Galaxy 


Visibility 


Beam size 


Beam size 


rms noise 


3cr column density 


3cr Hi mass 


ID 


weighting 


(arcsec) 


(kpc) 


(mJy beam~^) 


sensitivity (10^^ cm~^) 


sensitivity (10^ M©) 


UGC 2082 


uniform 


24.4 X 11.3 


1.7x0.8 


0.28 


5.4 


6.6 




robust-0 


29.1 X 12.6 


2.0 X 0.9 


0.21 


3.0 


4.9 




30'' taper 


40.9x33.1 


2.9 X 2.3 


0.20 


0.8 


4.7 


NGC 672 


uniform 


27.4 X 12.3 


1.0x0.5 


0.28 


4.4 


1.8 




robust-0 


30.6 X 13.5 


1.1 X 0.5 


0.21 


2.7 


1.4 




30" taper 


41.7x32.0 


1.5x1.2 


0.19 


0.8 


1.2 


NGC 925 


uniform 


22.2 X 12.6 


1.0x0.6 


0.24 


4.5 


2.3 




robust-0 


25.3 X 13.9 


1.1 X 0.6 


0.18 


2.7 


1.7 




30" taper 


37.9 x33.2 


1.7 X 1.5 


0.17 


0.7 


1.6 


NGC 4565 


uniform 


27.4 X 12.0 


1.4x0.6 


0.24 


3.9 


3.2 




robust-0 


31.9 X 13.8 


1.7x0.7 


0.17 


2.0 


2.3 




30" taper 


43.7 X 33.7 


2.3x1.8 


0.23 


0.8 


3.0 



The data cubes were used to generate images of the total H i 
line intensity (moment-0 images), and of the intensity-averaged 
line-of-sight velocity (first-moment images). These were pro- 
duced by smoothing the data cubes to a resolution of 90'', clip- 
ping the original cubes at the 3cr level of the smoothed cubes, 
and performing the moment calculations mentioned above. The 
miriad task moment was used for this process. In the case of 
the total intensity map, the primary beam correction was applied 
(for the WSRT 25 -m antennas, the primary beam shape takes the 
mathematical form cos^ (cry), where c is a constant with a value 
of 68 at L-band, v is the frequency in GHz, and r is the field 
radius in radians). Then, intensities were converted to column 
densities (A^m) assuming optically thin gas. The resulting maps 
are displayed and discussed in §|4] We note that typical veloc- 
ity profiles are more complicated than can be encapsulated in a 
single number such as the average velocity. Especially in edge- 
on galaxies, first-moment maps can be misleading, and therefore 
for our two highly inclined targets we instead show Renzograms 
(a single outer contour for each of a series of velocity channels, 
plotted on the same map). 

The primary goal of the survey is to characterize each galaxy 
in the sample by the properties of the extraplanar gas (if any) de- 
tected in the field, and to consider the full sample together to 
draw conclusions about the global characteristics of this phe- 
nomenon. A detailed picture of the individual galaxies will be 
one of the output products of this investigation. The results for 
such individual survey targets will be presented in a series of 
companion papers. The main engine for analyzing these galax- 
ies will be tilted-ring modeling software, which builds geomet- 
rical galaxy models using a series of idealized, roughly concen- 
tric, tilted rings. The tilted-ring modeling software is based in 
the GIPSY package. We will use a combination of GIPSY's 
GALMOD, and the T ilted Ring Fitting Code (TiRiFiC) which 
has been developed by Jozsa et al. ( 2007| ). 



4. Overview of Pilot Survey Observations 

In this Section, we present an overview of the data obtained dur- 
ing our Pilot Survey, which is made up of four targets observed 
during the first observing semester of the HALOGAS survey. 
The targets presented here are representative of a broad range of 
galaxy characteristics, and give a sense of the diversity in the full 
HALOGAS sample. For each galaxy, we produced a continuum 
map created using all 10 x 12hr of data, resulting in rms noise 
levels of about 30yuJybeam"^ All targets are detected in con- 
tinuum emission. We present maps of the total H i line intensity. 



comparisons with the optical image, and the Hi first-moment 
image or Renzogram. 

For each galaxy we provide some discussion of the target and 
the results gleaned from our HALOGAS Pilot Survey data prod- 
ucts. In the interest of brevity, we defer fuller detailed analyses 
of these and certain other individual HALOGAS targets (includ- 
ing NGC 4244; Zschaechner et al., in prep., UGC 4278; Serra et 
al., in prep., and NGC 949; Jozsa et al., in prep.) to forthcoming 
papers. Ultimately, the main goal of the HALOGAS project is to 
examine the properties of extraplanar H i across the full sample; 
those results will be presented following the completion of all 
survey observations. 

4.1. UGC 2082 

UGC 2082 is an edge-on, low surface brightness galaxy which 
appears to be isolated (hence the lack of environmental infor- 
mation in the UGC 2082 entry in Table [2|). It contains molecu- 
lar gas ( Matthews & Gao|[20nTt patthews et aLl[20Q5] ), but has 
a low star formation rate. The SFR estimate listed in Table [T] 
(0.023 Mq yr"^) is based on FIR luminosity, and by this estima- 
tor UGC 2082 has one of the lowest SFRs in our sample. We 
note, however, a confident detection of extended radio contin- 
uum emission from the inner disk. 

With the HALOGAS Hi observations, we detect a total Hi 
mass of Mhi = 2.54 x 10^ M©. The vertical distribution of Hi 
emission is somewhat extended, reaching up to a maximum ver- 
tical extent of about 5 kpc, fully one third of the radial extent of 
the Hi. The total Hi contours along the major axis clearly show 
a slight warp in the sky plane, beginning somewhat beyond the 
optical radius. A warp component along the line of sight may be 
responsible for some of the apparent thickening of the disk. 

Indeed, the appearance of the position-velocity (PV) dia- 
gram shown in Figure |2] is suggestive of a line-of-sight warp. 
In this diagram, the major axis is located at zero angular off'set. 
The characteristic triangular appearance, caused by emission at 
large vertical off'sets appearing at velocities close to systemic, 
can be caused by a halo with lagging rotationa l velocities or by 
a line-of-sight warp (e.g. 'Ooste rloo et al.||2Q0 7). In the case of 
UGC 2082, the visual appearance (in particular the narrowness 
of the velocity profiles at the largest angular off'sets) is more sim- 
ilar to galaxies determined to show evidence for a line-of-sight 
warp (e.g., |Rand & Benjamin]|20Q8[ [Gentile et al.||2QQ3| ) rather 
than those determined to host a lagging halo. The Renzogram 
shown in Figure [T] also shows characteristic signs of a line-of- 
sight warp: the contours of channels at intermediate velocities 
extend to larger radii with increasing distance from the mid- 



6 



G. Heald et aL: The WSRT HALOGAS Survey 



plane. That this effect is more pronounced on the approaching 
side, together with the uneven thickness of the H i distribution 
along the major axis, indicate that the line-of- sight warp is asym- 
metric. Geometrical modeling work is in progress to draw con- 
clusions about how much of this vertical extension can be due 
to a halo, ra ther than warping (o r flaring). Based on modeling 
as shown by Swaters et al. ( 1997 their Figure 4) and Qosterloo 
|et al. I ( [20071 their Figure 12), the diff'erent situations can be dis- 
tinguished by careful inspection of the individual channel maps, 
as well as PV diagrams such as the one shown in Figure [2] A 
detailed description of the modeling work is beyond the scope 
of this pilot survey description paper, and it will instead be pre- 
sented in a forthcoming publication. 

It is unlikely that large-scale morphological features such as 
warping or flaring can explain all of the high-latitude H i emis- 
sion, however. The "fingers" (filamentary features) of H i, ex- 
tending upward perpendicular to the major axis at the largest z 
on the northeast side, appear to be distinct enhancements of cold 
gas. One such feature is indicated with an arrow in Figure |2] At 
a projected distance of 9 kpc above the midplane, this cloud 
would have an Hi mass of a few 10^ M© and a broad velocity 
width of 60kms"^ The presence of features like this one sug- 
gests that at least some of the thick disk of UGC 2082 is made 
up of a filamentary H i halo. 

4.2. NGC 672 

NGC 672 is a moderately inclined spiral which has a neigh- 
boring irregular galaxy, IC 1727. In Figure [3] NGC 672 is 
the brighter galaxy centered in the frames, and IC 1727 is lo- 
cated to its southwest. At the adopted distance listed in Table 
[T] (7.6 Mpc), the linear projected separation between the two 
galaxies is 16.3 kpc. In the optical image, there is no clear sign 
of a strong d isturbance in NGC 672. Previous interferometric 
observations ( [Combes et al.|198Q ) have revealed a gravitational 
interaction with a substantial amount of gas in the interface re- 
gion. In our observations, we detect a total mass of 3.78x10^ M©, 
of which approximately 60% is associated with the primary, 
NGC 672. Both galaxies are clearly detected in continuum emis- 
sion. 

The H I distribution shown in Figure [3] reveals clearly the 
connection between the two galaxies. Each of the two galaxies 
appears to be relatively undisturbed within their optical extent, 
but the outskirts and interface region are more complicated. The 
pair seems to be in an early stage of interaction. A clear enhance- 
ment of Hi emission is located between the two galaxies, cen- 
tered at about aj2ooo.o = 0lM7'^4r, ^j2ooo.o = 27^22^^10^ To the 
northwest of the interaction region, a clear arc- shaped feature is 
seen extending away from the system. A faint diff'use extension 
is seen on the southwest side of IC 1727. A long, coherent arm- 
shaped feature starts just north of the center of NGC 672 and 
stretches away to the northeast. This could be a tidal or spiral 
arm, or may be a large Hi loop seen in projection. 

The first-moment image indicates that the rotation of 
NGC 672 is fairly well ordered within the optical body on the 
northeast side (farther from the interaction region), and signifi- 
cantly more disturbed toward the companion galaxy. The main 
body of IC 1727 appears to show well ordered rotation. The 
gas in the system that is located outside of the optical disks 
also has kinematics which overall blend smoothly into the first- 
moment images of the individual galaxies. It should be noted, 
however, that the velocity profiles are complicated in this region; 
this eff'ect is clearly visible in the PV diagram shown in Figure 
I?! (the emission near angular off'sets of 3' - 4'). There, the pro- 



files are significantly broadened, often have multiple peaks, and 
even show evidence for discrete emission at distinctly diff'erent 
velocities (e.g. the feature at angular off'set ~ 3\ and velocity 
420 - 460 km s"^). For that reason, the mean velocities in the 
interface region shown in Figure [3] should be regarded with sus- 
picion. The PV diagram also shows that the systemic velocity 
of IC 1727 is consistent with a continuation of the kinematics 



of NGC 672. As noted by Sancisi et al.[ ( g008) , in many minor 
mergers the kinematics of the infalling companion galaxy fol- 
lows that of the outer disk of the major galaxy. 

In order to gain a better understanding of this system, we 
intend to use Identikit (Barnes & Hibbard'2009 ), or similar soft- 
ware, to match the observations against computer models of 
galaxy interactions. This should give an indication of the three- 
dimensional geometry of the system, as well as the interaction 
history. There is clearly a large quantity of unsettled gas in the 
system, which is potentially material that will be accreted at a 
later stage. The ultimate fate of this gas is of obvious interest. 

4.3. NGC 925 



NGC 925 is a moderately inclined spiral. We detect bright dif- 
fuse continuum emission, with enhancements corresponding to 
the locations of the most prominent features in the optical spi- 
ral arms. NGC 925 has been observed previously in the Hi line, 
notably by Pisan o et al.[ ( [T998) and as part of THINGS (The 
Hi Nearby Galaxy Survey; [Walter et al.|2008] ). Neither of those 
studies detected H i at such large galactocentric radii as revealed 
in Figure [5] Note that the image shown by [Sancisi et al.|( [2Q08| ), 
based on earlier WSRT data, clearly illustrates the radially ex- 
tended and disturbed Hi distribution. With the HALOGAS ob- 
servations, we recover the same structure, but with higher sen- 
sitivity. The inner disk within the optical body has an unsmooth 
appearance, with clumps and holes. Many of the Hi overdensi- 
ties correspond to the locations of bright knots of optical light in 
the spiral arms. The appearance of the outer disk of NGC 925 
is that of a complicated stream of gas stripped during a minor 
merger. The twisting of the isovelocity contours in that region 
suggest outer disk warping or spiral structure. An enhancement 
in the Hi column density is seen at the southernmost edge of 
the outer disk, at coordinates aj2ooo.o = 02^26^^44^, ^j2ooo.o = 
33^25^^20^. Inspection of the optical DSS plates reveals that a 
faint fuzzy optical counterpart is very close to the location of 
this Hi enhancement, centered at about aj2ooo.o = 02^26^^53^, 
^j20oo.o = 33^25^^40^. Thus, the Hi streams may be the result 
of the disruption of a gas-rich dwarf companion, whose stellar 
counterpart is still barely visible in the survey plates. NGC 925 
is a prototypical minor merger, illustrating how gas is fed into 
large galaxies by the accretion of material from satellites. 

With the HALOGAS observations we detect 5.89x10^ M© of 
H I. A large amount of anomalous H i emission is clearly visible 
in this target. While the streams in the outer disk generally follow 
the overall rotation of the inner disk (see the first-moment image 
in Figure [5]), there is also gas in the in ner parts of the galaxy 
which either appear as "beard" emission ( [Fraternali et al.|2002 ), 
or sometimes even at distinct velocities relative to the bulk of the 
galaxy rotation. These features cannot be seen directly in the mo- 
ment maps shown in Figure [5] The beard gas is however readily 
visible in the PV diagram shown in Figure |6| There, faint emis- 
sion is seen to extend preferentially toward the systemic velocity 
(to lower velocity at positive angular off'set, and higher velocity 
at negative angular off'set). This may well be the sign of a slowly 
rotating halo seen in projection against the disk. Interpreted in 
this way, we estimate the total H i mass of the lagging compo- 
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Fig. 1. Overview of the HALOGAS observations of UGC 2082. The left panel shows the Hi total intensity overlaid on the DSS 
7?-band image. The Hi contours originate from the 3C-tapered image, begin at Nm = 1.0 x 10^^ cm"^ and increase by powers of 
two. The straight line shows the orientation of the PV slice shown in Figure [2] The right panel shows an overlay of several channels 
in the lowest resolution data cube, all at a level of 0.9 mJy beam" ^ 3.75cr). The contours are separated by 12.4kms"^ begin at 
593 km s~^ (dark blue) and range upward to 815 km s"^ (dark red). Both panels show the same area of the sky. The beam size of the 
H I data is shown in the lower left corners of the left panel. 












— - — :=> 



Angular offset [arcmin] 

Fig. 2. Position- velocity (PV) diagram for UGC 2082, plot- 
ted using the 30'' -resolution cube. Contour levels begin at 
300yuJybeam"^ 1.5(t; see Table [3]) and increase by powers 
of two. The dashed contours are at -300yuJy beam"^ The line 
in Figure [T] indicates the orientation of the PV slice shown here. 
The northeast end of the slice is to the left. At the distance listed 
in Table[T] V corresponds to a linear scale of 3. 1 1 kpc. The arrow 
indicates an extraplanar feature described in the text. 



nent to be a few times 10^ M©, at a velocity difference of about 
20 - 50kms"^ with respect to the disk. We note that the SFR 
of NGC 925 is similar to that of NGC 2403, which is the pro- 
totypical galaxy demonstrating beard emission, and the latter is 
somewhat more massive. Detailed modeling of NGC 925 will al- 
low us to separate the regularly rotating disk emission from the 
anomalous gas, determine its nature, and make a more detailed 
comparison with NGC 2403. 



A small companion (outside of the frame in Figure [5]) is 
found to the north of NGC 925 at about Qfj2ooo.o = 02^27'^20^ 
^j20oo.o = 33^57^^30^ Its Hi mass is 1.27 x 10^ M©. This object, 
too, appears to have a faint optical counterpart in the DSS plates. 
It is located at a projected distance of 60 kpc from the center of 
NGC 925 and does not show any obvious signs of interaction. 



4.4. HGC 4565 

NGC 4565 is a large, edge-on spiral galaxy which is similar to 
the Milky Way. The new HALO GAS observations ar e combined 
with archival data presented by Dahlem et al. ([2005 ). The other 
two galaxies in the frame are NGC 4562, a small companion 
classified in the RC3 as type SB(s)dm and visible to the south- 
west of NGC 4565, and IC 3571, the small galaxy located very 
close to NGC 4565 on its north side. In the PV diagram shown in 
Figure [8j the vertically extended emission in NGC 4565 is seen 
to have a smooth velocity gradient directly toward the systemic 
velocity of IC 3571. 

The distribution of H i emission is shown in Figure [7] The 
total H I mass detected in NGC 4565 itself is 7.43 x 10^ M©. The 
galaxy is obviously warped: this can be seen in the total H i map 
(particularly obvious on the north end of the disk) and also in the 
first-moment image contours. NGC 4562 has a regularly rotating 
Hi disk with a total Hi mass of 1.82 x 10^ M©. IC 3571 is also 
detected in Hi (Mm = 4.18 x 10^ M©), and the Hi data show 
clear signs of a minor interaction between it and NGC 4565. 



NGC 4565 has a moderately high star formation rate when 
compared to the rest of the HALOGAS sample. The most similar 
HALOGAS sample galaxy (with existing observations to date) is 
NGC 891. There are, however, several key diff'erences between 
these two galaxies. First, the SFR of NGC 4565 is only about 
one quarter of that of NGC 891. Both galaxies have high rota- 
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Fig. 3. Overview of the HALOGAS observations of NGC 672. The left panel is as in Figure T] The contours begin at 
2.5 X 10^^ cm"^ and increase by powers of two. The right panel shows the first-moment image derived from the lowest-resolution 
H I cube. The grayscale ranges from 245 km s"^ (black) to 570 km s"^ (white). Contours are white on the approaching side and black 
on the receding side, and run from 255 km s"^ to 600 km s"^ in steps of 15kms"^ The white contour closest to the center of the 
primary galaxy, corresponding to 425 kms"\ represents the systemic velocity of NGC 672. 




Angular offset [arcmin] 



Fig. 4. As in Figureji) but for NGC 672. The line in Figure 3]in- 
dicates the orientation of the PV slice shown here. The northeast 
end of the slice is to the left. At the distance listed in Table [T] V 
corresponds to a linear scale of 2.18 kpc. 



tion speeds (but that of NGC 4565 is 15% higher - it has the 
highest rotation speed of all the HALOGAS sample galaxies). 
The distances to NGC 4565 and NGC 891, and their optical di- 
ameters, indicate that NGC 4565 is larger in linear dimensions. 

It is apparent from Figure [7] that the H i distribution in this 
galaxy is far less ver tically extended than that of NGC 891 (com- 
pare with Figure 1 of |Oosterloo et al.|2007| ), even though our new 
observations have comparable sensitivity to the NGC 891 work. 
Although the disk of NGC 4565 is clearly somewhat warped 
(and importantly, also partially along the line of sight; ^Rupen, 
|1991| ), we can still derive a reasonable upper limit to t he verti- 



cal extent, of approximately 5 kpc (using RA = 315°; |Dahlem 
[eTa l. 2005). Meanwhile, the Hi halo of NGC 891 extends up 
to more than 10 kpc - a distinct diff'erence. Previous work has 
also determined that there is a lack of ionized gas in the halo of 



NGC 4565 (see |Rand etal.|1992 i A final diff'erence between the 
two galaxies is that NGC 891 has a extremely bright and verti- 
cally extended radio continuum halo (e.g., pahlem et al.,1994j ). 
We find that NGC 4565 also has a vertically extended, and to 
our knowledge previously undetected at these frequencies, radio 
continuum halo (see Figure |9]). However, it does not seem to be 
as extended (the continuum halo in NGC 4565 reaches up to 2-3 
kpc, compared to 10 kpc in NGC 891), nor is it as bright, as the 
halo in NGC 891. Studies of the continuu m emission (and its po- 
larization) in NGC 4565 are d escribed by Hummel et al.| ( |1984| ) 
and Suku mar & Allen] ( |1 99 1| ). Further discussion regarding the 
diff'erences between NGC 4565 and NGC 891 is provided in §[5] 



5. Discussion 

The HALOGAS pilot survey, while still limited in the number 
of targets, nevertheless covers a substantial range of galactic 
properties. From the perspective of environmental properties, we 
have an apparently isolated galaxy (UGC 2082), two targets that 
have small, gas-rich companions (NGC 925 and NGC 4565), 
and one galaxy in an interacting pair (NGC 672). We also sam- 
ple about two orders of magnitude in star formation rate, and 
nearly a factor of three in rotation speed, with this pilot survey. 
While it is clear that none of the pilot survey galaxies have an 
Hi halo of the scale seen in NGC 891, we find that all of the 
pilot survey targets show varying indications of cold gas in their 
halos. 

The lack of a prominent halo in UGC 2082 may not be a 
surprising result: the star formation rate is low, and the galaxy 
is isolated. Apart from the possibility of significant primordial 
accretion, a ready source for creating a gaseous halo does not 
seem to be present in this target. Nonetheless, we find several 
distinct extraplanar features (these are noted in § 4.1 and are vis- 
ible in Figures [T] and |2]). These clouds may be a sign of cold gas 
accretion. 
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Fig. 5. Overview of the HALOGAS observations of NGC 925. The panels are as in Figure [s] The contours in the left panel begin at 
Nhi = 2.0 X 10^^ cm"^ and increase by powers of two. The grayscale in the right panel ranges from 440 km s"^ (black) to 675 km s"^ 



(white). Contours are white on the approaching side and black on the receding side, and run from 450 km s ^ to 660 km s ^ in steps 
of 15 km s"^ The white contour closest to the center of the galaxy, corresponding to 555 km s"\ represents the systemic velocity of 
NGC 925. 




Angular offset [arcmin] 



Fig. 6. As in Figure (i) but for NGC 925. The line in Figure [i] 
indicates the orientation of the PV slice shown here. The east 
end of the slice is to the left. The lowest contour level (solid and 
dashed) is ±250 yuJy beam" ^ At the distance listed in Table[l] V 
corresponds to a linear scale of 2.73 kpc. 



Two of the pilot galaxies that are located in denser environ- 
ments have a much less regular H i structure. Clearly, the case 
of NGC 672 is a complicated one - the strong interaction with 
its companion IC 1727 has spread Hi emission throughout the 
system. This is a dramatic case of interaction-induced accretion. 
Another interacting system but at a later stage, NGC 925 seems 
to have interacted with a gas-rich companion sometime in its 
recent past, resulting in stripped material surrounding the main 
disk. Whether the culprit is the faint blob at the southern ed ge o f 
the galaxy, the object found to the north (as described in § 4.3 ), 
or some other entity, is not yet clear. NGC 925 also seems to 
possess H I at lagging velocities with respect to the main disk. 



The anomalous component in NGC 925 has a similar mass and 
kinematic lag when compared to the same feature in NGC 2403. 



Another target to compare with existing observations is the 
large edge-on NGC 4565, which does not show evidence of a 
prominent halo component. Can this lack be explained? In com- 
parison with NGC 891, the differences in the global parameters 
of the two galaxies all tend to be in the sense in which one would 
indeed expect a larger halo in NGC 891, if gaseous halos are in- 
duced primarily by star formation-driven galactic fountain flows. 
Under that assumption, the most important properties of a galaxy 
for determining whether a halo is formed should be (1) the SF 
energy injection rate (per unit disk area; see e.g. Dahl em et al.| 
|1995HRossa & Dettmar 2003 ), and (2) the gravitational poten- 
tial (traced by v.^^^). WithNGC 891 's larger SFR, smaller disk 
area (at the distances listed in Table [TJ, and lower rotational ve- 
locity, it would seem that it should indeed be more efficient in 
maintaining a gaseous halo. Per unit disk area, more star forma- 
tion energy is available, and gravity is less eff'ective at prevent- 
ing material from moving away from the disk midplane. That 
NGC 925 and NGC 2403 have similar anomalous gas features, 
as well as similar SFR and mass, would also be consistent with 
such a picture of gaseous halo formation. 



Do the other HALOGAS targets follow the same trend with 
the global galaxy parameters? Though it is intriguing, the num- 
ber of galaxies which have been observed with enough sensitiv- 
ity to distinguish a population of halo gas remains rather low, so 
that general statements are still uncertain. This initial result from 
our pilot survey underlines the need for a suitably large survey, 
in order to trace a range of variables in the population of nearby 
galaxies. The full HALOGAS survey will provide such a large 
sample. 
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Fig. 7. Overview of the HALOGAS observations of NGC 4565. The panels are as in Figure[l] The contours in the left panel begin at 
Nhi = 2.0 X 10^^ cm"^ and increase by powers of two. The contours in the right panel are all at 1.0 mJy beam" \ and are shown for 
a number of channels separated by 12.4 km s"\ starting at 945 km s"^ (dark blue) and ranging upwards to 1489 km s"^ (dark red). 




4 2 -2 -4 

Angular offset [arcmin] 



Fig. 8. As in Figure [2) but for NGC 4565. The line in Figure [t] 
indicates the orientation of the PV slice shown here. The north- 
east end of the slice is to the left. The lowest contour level (solid 
and dashed) is ± 350 yuJy beam" ^ At the distance listed in Table 
[T] r corresponds to a linear scale of 2.82 kpc. 

6. Multi-wavelength supplements 

For a full census of the contents of the survey galaxies, our 
H I observations must be supplemented with high-quality, sen- 
sitive observations tracing the distribution and kinematics of 
stars, dust, and other gas phases. With that in mind, the WSRT 
HALOGAS program is being supplemented in several ways. 

Since minor mergers are thought to be one of the major chan- 
nels for gas accretion, stellar- and gas accretion are profoundly 
linked and cannot be studied independently. We therefore ob- 
tained observing time with the Isaac Newton Telescope (INT; 
using the Wide Field Camera) in order to complement the deep 
HALOGAS H i observations with a photometric survey suited 
to detect the stellar remnants of former accretion events. The 
INT observations for a subset of the HALOGAS sample have 
been obtained, and the data will be presented in a forthcom- 




RA 02000) 



Fig. 9. Radio continuum image of NGC 4565. The contour levels 
begin at 54yuJy beam"^ 3cr) and increase by powers of two. 



ing companion paper. This optical observing program is named 
HALOSTARS. 

The 3.5m ARC telescope at Apache Point Observatory is be- 
ing used for kinematic observations of the ionized gas in a sam- 
ple which overlaps the edge-on subsample of HALOGAS, us- 
ing a novel multi-slit technique that allows for simultaneous Ha 
spectra of 1 1 parallel long slits, spaced 22'' apart over a 4.5' x4. 5' 
field of view. Together this set of long slits can be used to form 
a velocity field. So far, data have been obtained for NGC 891, 
NGC 3044 (not a HALOGAS target), NGC 4244, NGC 4565, 
and NGC 463 1 . The observations will continue through the next 
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several observing seasons. These optical observations will allow 
us to test whether both the ionized and neutral phases share the 
same origin a nd kine matics (as has been previously investigated 
in NGC 891; |Heald e t al. 2006 ). Moreover, the ionized compo- 
nent is more tightly coupled to the radial extent of the star form- 
ing disk than is the neutral, likely resulting in little to no contam- 
ination from warped outer disk gas complicating the analysis of 
extraplanar emission. This provides an important check on the 
H I modeling results. 

We have received 102 ksec on the Galaxy Evolution Explorer 
(GALEX; [Martin et al.|2005| ) satellite to observe 8 galaxies, of 
which 6 are drawn from the HALOGAS sample and 2 from 
a complementary Very Large Array (VLA) Hi program. In 
addition, 3 HALOGAS targets are already being observed as 
part of the GALEX Deep Galaxy Survey (DGS) survey, while 
NGC 2403 and NGC 891 also have existing deep archival 
GALEX exposures from other programs. We were granted ob- 
servations in both GALEX bands, but unfortunately the FUV 
channel onboard GALEX is now unavailable, so that we will 
only obtain NUV images. The purposes of the GALEX obser- 
vations are (/) to search for scattered light from dust co-located 
with the gaseous halo in edge-on systems, which would help elu- 
cidate the origin of this material; and (//) to perform a highly 
sensitive search for star formation in the outer disks and halo re- 
gions in all of the galaxies, allowing us to study the correlation 
between gas accretion and outer disk star formation. 

Finally, we also plan on obtaining wide-field UBRHck imag- 
ing in support of the GALEX and WSRT observations to search 
for outer disk and halo star formation, and for companions to de- 
termine to what extent accretion of gas may be associated with 
gas rich faint dwarfs. 

7. Data release and availability 

The calibrated survey data will be made publicly avail- 
able in stages as the observations, data reduction, and 
analysis progress. The survey webpage is located at 
http://www.astron.nl/halogas and this is the primary 
location at which reduced data cubes, moment maps, and other 
data products can be found. Data releases will coincide with the 
publication of analysis papers. Pre-release data products may be 
obtained by contacting the HALOGAS team. 
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